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ABSTRAGT 

Using a sample of dwarf galaxies observed using the VIMOS IFU on the VLT, we investigate the 
mass-metallicity relation (MZR) as a function of star formation rate (FMRsfr) as well as Hl-gas 
mass (FMRhi). We combine our IFU data with a subsample of galaxies from the ALFALFA HI 
survey crossmatched to the Sloan Digital Sky Survey to study the FMRsfr and FMRhi across the 
stellar mass range 10^’^ to 10^’^ Mq, with metallicities as low as 12+log(0/H) = 7.67. We find the 
Icr mean scatter in the MZR to be 0.05 dex. The la mean scatter in the FMRsfr (0.02 dex) is 
significantly lower than that of the MZR. The FMRsfr is not consistent between the IFU observed 
galaxies and the ALFALFA/SDSS galaxies for SFRs lower than 10“^‘^ Mq yr“^, however this could 
be the result of limitations of our measurements in that regime. The lowest mean scatter (0.01 dex) 
is found in the FMRhi- We also find that the FMRhi is consistent between the IFU observed dwarf 
galaxies and the ALFALFA/SDSS crossmatched sample. We introduce the fundamental metallicity 
luminosity counterpart to the FMR, again characterized in terms of SFR (FMLsfr) and Hl-gas mass 
(FMLhi). We find that the FMLhi relation is consistent between the IFU observed dwarf galaxy 
sample and the larger ALFALFA/SDSS sample. However the la scatter for the FMLhi relation is not 
improved over the FMRhi scenario. This leads us to conclude that the FMRhi is the best candidate 
for a physically motivated fundamental metallicity relation. 

Key words: galaxies: abundances - galaxies: formation - galaxies: starburst - galaxies: dwarf - 
galaxies: evolution 


1. INTRODUCTION 

The mass-metallicity relation (MZR) is an impor¬ 
tant tool for understanding the underlying processes of 
galaxy formation. A galaxy’s metal content provides in¬ 
sights into the history of star formation within a galaxy. 
Lequeux et al. (1979) and Kinman & Davidson (1981) 
first reported on the correlation between galaxy stellar- 
mass and gas-phase metallicity. Later studies confirmed 
the relation using the Sloan Digital Sky Survey (SDSS) 
for large number statistics (Tremonti et al. 2004; Kewley 
& Ellison 2008; Michel-Dansac et al. 2008; Salim et al. 
2014). Other studies have focused on the redshift evolu¬ 
tion of the MZR (Savaglio et al. 2005; Erb et al. 2006; 
Cowie & Barger 2008; Liu et al. 2008; Maiolino et al. 
2008; Panter et al. 2008; Rodrigues et al. 2008; Hayashi 
et al. 2009; Mannucci et al. 2009; Perez-Montero et al. 
2009; Cresci et al. 2012). 

The luminosity-metallicity relation (LZR) has also 
been proposed as a relation of merit (Skillman et al. 1989; 
Garnett 2002; Perez-Gonzalez et al. 2003; Pilyugin et al. 
2004; Lee et al. 2006). The advantage of the LZR is that 
reliable luminosity measurements are generally easier to 
obtain than reliable stellar mass estimations. The dis- 
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advantage of the LZR is that it is dependent upon the 
bandpass being used, and could be affected by dust ex¬ 
tinction (Salzer et al. 2005). The LZR is observed to 
hold over a range of 10 magnitudes in galaxy optical lu¬ 
minosity (Zaritsky et al. 1994; Tremonti et al. 2004; Lee 
et al. 2006). The MZR and the LZR reveal that some 
physical mechanism must be driving a correlation be¬ 
tween stars, metals, and gas flows across a vast range 
of scales. However, the exact physical mechanism(s) re¬ 
sponsible for forming this relation are still uncertain. 

Recent studies have investigated the mass-metallicity 
relation with a particular focus on low stellar-mass or 
low-luminosity systems (Berg et al. 2012; Nicholls et al. 
2014; Haurberg et al. 2015; McQuinn et al. 2015; James 
et al. 2015), which is the population of galaxies of inter¬ 
est to this work. The MZR is known to exhibit scatter 
greater than would be expected by the uncertainties on 
the individual data points (Tremonti et al. 2004; Lara- 
Lopez et al. 2010; Mannucci et al. 2010). 

To explain this apparently enhanced scatter, it has 
been suggested that the MZR has a dependence on Star 
Eormation Rate (SER; Lara-Lopez et al. 2010; Mannucci 
et al. 2010; Andrews & Martini 2013), which could alter¬ 
natively be explained as a dependence on Hl-gas content 
(Bothwell et al. 2013; Lara-Lopez et al. 2013). This 3 
dimensional extension of the mass-metallicity relation is 
known as the fundamental metallicity relation (EMR). 
The EMR has been shown to be consistent, with little 
evolution, up to a redshift of 2.5 (Lara-Lopez et al. 2010; 
Mannucci et al. 2010). It is currently uncertain as to 
which observable (HI mass or SER) better explains the 
scatter in the MZR, or if some other unknown parameter 
is responsible. 

Several hypotheses have been proposed for the mech- 
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anisms that guide the FMR, such as pristine gas inflows 
diluting metal abundances (Finlator & Dave 2008; Dave 
et ah 2010). Inflows of pristine gas would cause a galaxy 
to exhibit lower gas-phase metallicity than would be pre¬ 
dicted by the FMR. Recent evidence suggests that infall 
of pristine gas is the more dominant component at high 
redshift (Agertz et al. 2009; Bournaud & Elmegreen 2009; 
Brooks et al. 2009; Dekel et al. 2009) 

As argued in Mannucci et al. (2010) infalling pristine 
gas must eventually ignite star formation, which will 
create strong supernova winds that will expel metals. 
Gas outflows expelling metal-rich gas into the interstellar 
medium provide an alternative explanation for the ob¬ 
served gas-phase metallicity being lower than the MZR 
prediction (Edmunds 1990; Lehnert & Heckman 1996; 
Garnett 2002; Tremonti et al. 2004; Kobayashi et al. 
2007; Scannapieco et al. 2008; Spitoni et al. 2008). 

Studies on the EMR thus far, either as a function of 
SER (EMRsfr) ot hi mass (EMRhi), have focused on 
the SDSS sample of spectroscopically observed galaxies 
(Lara-Lopez et al. 2010; Mannucci et al. 2010; Bothwell 
et al. 2013). As such these galaxies are biased towards 
the intermediate to high stellar-mass galaxies that dom¬ 
inate the SDSS spectroscopic sample. However it is im¬ 
portant to understand the low-mass galaxy population as 
well and to test whether or not they fall onto the same 
EMR. 

Dwarf galaxies are the most common type of galaxies 
in the Universe, comprising ^ 85% of the objects in the 
volume within lOMpc (Karachentsev et al. 2004). Even 
though they are not as widely studied as their more mas¬ 
sive counterparts, they may hold the answers to several 
fundamental questions concerning the processes of galaxy 
formation and evolution. 

The data presented here consists of integral field unit 
(lEU) observations of a sample of 11 dwarf galaxies se¬ 
lected from the ALEALEA survey. ALEALEA is a blind 
HI survey of the local universe {cz <180,000 km s“^; 
Haynes et al. 2011). The lEU observations were taken 
using the VIMOS lEU spectrograph on the Very Large 
Telescope (VLT). 

To compliment the lEU observed dwarf galaxy dataset, 
we utilize a larger ALEALEA/SDSS cross-matched sam¬ 
ple of spectroscopically observed galaxies. In addition we 
also include samples from the literature of low-mass/low- 
luminosity galaxies such as the Berg et al. (2012), long- 
slit spectroscopic observations of low-luminosity galaxies, 
the SHIELD survey of ALEALEA dwarf galaxies (Hau- 
rberg et al. 2015; McQuinn et al. 2015) and the James 
et al. (2015) long-slit morphologically selected survey of 
dwarf irregular galaxies. 

Throughout this paper, we assume a Hubble constant 
of Ho = 67 km s“^ Mpc“^, Qm = 0.32, and Qa = 0.68. 


2. SAMPLE 

The lEU observed dwarf galaxy sample includes nearby 
(D < 20Mpc), low HI mass (Mhi < 10^‘^), low surface 
brightness dwarf galaxies, selected to study the mass- 
metallicity relation on the low stellar-mass end. Due to 
the nature of observing nearby dwarf galaxies, lEU spec¬ 
troscopy is vital to capture the properties of the entire 


galaxy. As can be seen in Eigure 1, for the galaxies in 
our sample that have Sloan Digital Sky Survey (SDSS) 
fiber observations, the single SDSS fiber only captures 
a small portion of the galaxy’s flux, and the fiber may 
miss bright star forming regions (e.g. Richards et al. 
2014), this would cause star formation rate estimations 
of nearby dwarf galaxies such as these to be inaccurate. 

lEU observations also have the advantage of not re¬ 
quiring additional photometric observations to identify 
bright HII regions to target. Performing such observa¬ 
tions using long-slit spectroscopy would require a-priori 
knowledge of the positions of HII regions which can be 
obtained from narrow-band Ha images, introducing an 
optical bias in the otherwise Hl-selected sample. The 
most metal poor galaxies are also likely to have the 
faintest HII regions, making them barely detectable in 
narrow-band images. Integrating over a larger surface 
area allows us to target fainter galaxies. Because the 
dwarf galaxies that we are trying to observe are so faint, 
only 6 galaxies in this sample have reliable SDSS spec¬ 
troscopy. A purely SDSS based study, or a long-slit spec¬ 
troscopy study are both likely to preferentially exclude 
the faintest dwarf galaxies. 

2.1. Spectroscopic Observations 

Spectroscopic data of 28 dwarf galaxies were initially 
taken using the VIMOS (Le Eevre et al. 2003) lEU spec¬ 
trograph on the Very Large Telescope (VLT) located at 
Paranal Observatory. Of the 28 observed galaxies, 11 
have Ha emission lines greater than our amplitude over 
noise (AoN) cut of 5 for more than 20 spaxels and there¬ 
fore are included in this study. Integration times for the 
remaining 17 galaxies were insufficient to reach our tar¬ 
get depth. The data was obtained starting on April 11, 
2008 and ending on May 19, 2010 under program IDs 
081.B-0649(A) and 083.B-0662(A). Data was obtained 
using the VIMOS Low Resolution Blue Grism which has 
a wavelength range of 4000-6700 A and a spectral reso¬ 
lution of 5.3 A pixeU^ (R ^ 1000). 

Using the LR Blue grism provides the full 54"x54" field 
of view possible with VIMOS which allows us to obtain 
in a single pointing spectra across the complete stellar 
disk of each galaxy. Each object was observed using a 3 
dither pattern, with each dither being integrated for 20 
minutes. Average seeing across all observations is 1.05" 
EWHM. 

The LR blue grism provides both a wider spectral 
range, and a wider field of view when compared to the 
HR blue grism. The LR grism wavelength range allows 
for simultaneous observations of the emission lines from 
HP (A = 4861) to [Nil] (A = 6583). The major drawback 
to using the low-resolution spectra (^5.3 A pixeU^) is 
that we are unable to measure gas kinematics from the 
emission lines and the instrumental dispersion causes the 
Ha and [Nil] A A = 6549,6583 A emission lines to blend 
together. 

The spectroscopic data obtained with the VIMOS lEU 
is reduced from its raw form using the Reflex environ¬ 
ment for ESO pipelines (Ereudling et al. 2013). The 
standard VIMOS template is used within the Reflex en¬ 
vironment to produce the master bias and calibration 
frames containing the fiber traces and wavelength solu¬ 
tion. Many of the raw data and calibration frames con- 
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Fig. 1.— SDSS r-band images of the 11 galaxies presented in this study, demonstrating the need for IFU spectroscopy. The red circle 
indicates the location of the 3" SDSS fiber which would be used to obtain spectral information without IFU observations. The SDSS fiber 
often misses a considerable amount of light from the clumpy star-forming regions of these nearby dwarf galaxies. Shown in blue is the 
Petrosian radius as determined by the SDSS automated data reduction procedure. The SDSS photometry pipeline visibly overestimates 
or underestimates the flux in galaxies from this sample. Each thumbnail image is the size of the full VIMOS low resolution multiplexed 
field-of-view (54"x54" ). The green boxes indicate the subset of this field-of-view that our IFU analysis focuses on. 


tain a bright artifact across the surface of the chip, iden¬ 
tified to be an internal reflection within the instrument. 
This contamination interferes with the wavelength cali¬ 
bration routine within Reflex because it is often misiden- 
tified as a skyline, causing the spectrum to be shifted 
incorrectly. In order to compensate, we disable the sky¬ 
line shift in the calibration steps. The wavelength so¬ 
lution provided by the ESO pipeline proved to be ac¬ 
curate without using the skyline shift. We also use the 
flux standardization routine within the Reflex VIMOS 
pipeline. We then apply the calibrations frames to the 
science frames to produce the fully reduced Row-Stacked 
Spectra (RSS). 

The final output of the Reflex pipeline is four quad¬ 
rants per observation. We input these individual RSS 
quadrants into routines written in IDL and Python^. 
These routines fit Gaussian curves to the skylines and 
subtract the sky background while retaining the flux in¬ 
formation from each skyline to normalize the transmis¬ 
sion of each spaxel across the entire field of view. After 
the normalization and sky subtraction steps are com¬ 
pleted, the two dimensional RSS are converted into three 


dimensional data cubes. These data cubes consist of two 
spatial dimensions with the wavelength axis being the 
third dimension. Further details of the process these 
scripts follow can be found in Jimmy et al. (2013). 

Once the data cubes have been built, the individual 
dithers (typically 3 per galaxy) are stacked using a 5cr 
clipped mean. We use the AoN to select for spaxels con¬ 
taining sufficient emission line flux for oxygen abundance 
estimations. AoN is defined as the amplitude of the emis¬ 
sion line divided by the noise after the linear offset is 
subtracted. An AoN threshold of 5 is used to identify 
spaxels which contain emission lines for gas-phase metal- 
licity analysis. Only spaxels which pass the AoN cut in 
all of the three following lines are included in the inte¬ 
grated spectrum for oxygen abundance estimations: H/3 
A = 4861, [OIII], A =5007, Ra X = 6563. We do not 
require [Nil] A = 6583 to be separately detected because 
it is blended with Ra (Figure 2). 

All spaxels which pass our AoN cut are then fed into 
Python based Gaussian fitting routines to measure the 

^ Available publically: http://jimmy.pink/^^code 
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Fig. 2.— Example spectra showing the blending within the [Nil] 
and Ha emission lines due to the instrumental resolution (7.5 A). 
Our Python based Gaussian fitting routines are able to deconvolve 
the three emission lines into their constituent parts, as can be seen 
by the three Gaussian curves plotted in red, green, and blue. The 
flux in the blue [Nil] 6549 A emission line is negligible. Details of 
our ability to recover the three emission line fluxes can be found in 
Appendix A. Our oxygen abundance estimations rely upon mea¬ 
surements of the Ha and [Nil] emission lines. 

emission line flux ratios for gas-phase metallicity estima¬ 
tions. To obtain the integrated spectrum used to derive 
oxygen abundances, we sum the spectra from each spaxel 
that passes the AoN cut to produce a single spectrum 
that we perform this analysis upon. We chose to write 
our own routines instead of using a publicly available 
program such as GANDALF (Sarzi et al. 2006). Because 
GANDALF fits the stellar continuum and emission lines 
simultaneously, GANDALF is unable to fit robustly the 
weak stellar continuum in the dwarf galaxies. 

Using an AoN cut is necessary because a significant 
fraction of the spaxels in the data cubes are excessively 
noisy, and including them in the integrated spectrum 
would significantly reduce the signal-to-noise ratio. Un¬ 
fortunately the AoN cut causes us to discard a number 
of spaxels containing Ha flux, resulting in an underes¬ 
timate of the SFR. In order to ensure that we capture 
all of the Ha flux, we use the segmentation images de¬ 
scribed in Section 2.2 to select spaxels belonging to a 
galaxy. We rescale and rotate the segmentation map to 
match the VIMOS data cube using the IDL routine HAS- 
TROM.PRO^. The integrated spectra obtained using the 
segmentation map selection is used to measure Ha flux 
for the purposes of determining the SFR for each galaxy. 

2.1.1. Emission Line Fitting 

The basic procedure that our Python based Gaussian 
fitting routines follows is as such: First a linear fit to 
the spectrum continuum is performed using the opti- 
mize_curve routine to fit a line using a least squares 
optimization. The fitting window around each emission 
line is approximately 400 A (^75 pixels) with nearby 
emission, sky lines, and internal reflections regions being 
masked out. After the continuum has been subtracted, 

^ Available as part of the IDL Astronomy User’s Library 



the [Nil] and Ha emission lines are fit using gaussfitter^ 
to fit 3 Gaussian curves simultaneously. This allows us 
to separate the individual components of the 3 blended 
emission lines. 

The output measurements of gaussfitter are not heavily 
dependent upon initial input parameters assumed, there¬ 
fore we provide it an initial input intensity of 0.5 erg 
cm“^s“^, an initial input peak position of the redshift 
corrected expected peak for each line, and an initial in¬ 
put dispersion of 7.5 A, equal to the instrumental disper¬ 
sion of the observations. No line flux ratios are assumed 
within our fits. The gaussfitter bounds are set up such 
that each Gaussian fit is limited to a maximum 10% vari¬ 
ation from its expected peak position and FWHM. Be¬ 
cause the instrumental dispersion dominates the width of 
the emission lines, we do not expect line widths to vary 
widely across spaxels. 

We do not include a fit to the Balmer absorption in our 
procedures. Our observations are insufficient to reliably 
estimate the stellar continuum for accurate absorption 
measurements. However it is likely that Balmer absorp¬ 
tion effects are negligible and well within our uncertain¬ 
ties (Rosa-Gonzalez et al. 2002). Further details about 
the emission line flux fitting code and our ability to re¬ 
cover the blended emission lines can be found in Ap¬ 
pendix A. 

Uncertainties on our flux measurements are estimated 
using a Monte Garlo Markov chain (MCMC) technique 
with 1000 iterations. For each iteration, random noise is 
added to the spectrum on the order of the residuals to the 
linear fit. Then the emission line fluxes are remeasured. 
The la width of the Gaussian distribution of flux values 
provided by 1000 iterations provides the uncertainties in 
our line fluxes. These uncertainties are then propagated 
using standard Taylor series error propagation methods 
in order to determine the uncertainties in our derived 
quantities. 

We perform the linear and Gaussian fits on the blue 
and the red end of the spectrum separately, with the 
blue end containing the H^ and [OIII] lines, and the red 
end containing the Ha and [Nil] lines. The results of our 
fitting routine on the red end of an integrated spectrum 
can be seen in Figure 2. 

2.2. Photometric Observations 

To obtain stellar mass estimations for our galaxies, we 
utilize the stellar mass estimation procedure outlined in 
West et al. (2010). Doing so requires accurate photo¬ 
metric observations from the Sloan Digital Sky Survey 
(SDSS; Eisenstein et al. 2011) Data Release 12 (Alam 
et al. 2015) which is the final data release of SDSS- 
III. The automated SDSS photometric pipeline (Lupton 
et al. 2002) is optimized for small, well resolved objects 
and is known to be unreliable for clumpy and diffuse ob¬ 
jects such as dwarf galaxies, mostly due to the SDSS algo¬ 
rithm’s aggressive parent-child splitting routines. As can 
be seen in Figure 1, Petrosian flux measurements could 
possibly be an inaccurate representation of the flux of 
our galaxies. West et al. (2010) found that roughly 25% 
of SDSS galaxies have less than 90% of their flux con¬ 
tained in the brightest child. Therefore, to obtain more 

^ written by Adam Ginsburg 

http: / / www.adamgginsburg.com/pygaussflt.htm 
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accurate stellar mass estimations, we have chosen to run 
Source Extractor (Bertin & Arnouts 1996) manually on 
each of the SDSS bandpass image frames for each of the 
11 galaxies within our IFU sample. 

To ensure that the SDSS images are aligned prop¬ 
erly, we first use HASTROM.PRO to align each indi¬ 
vidual bandpass fits image to the r-band image coor¬ 
dinates. We stack all 5 photometric bands to create 
a higher signal/noise detection image. We run Source 
Extractor in dual image mode with the detection im¬ 
age and each of the individual frames. The default val¬ 
ues and keyword searches are used as inputs to Source 
Extractor except the following: MAG_ZEROPOINT 
= 22.5, BACK_SIZE = 256, and PIXEL_SCALE = 
0.396. Additionally for every galaxy we begin with 
the following parameters: DEBLEND_NTHRESH, DE- 
TECT_THRESH, ANALYSIS_THRESH set to 1.0, and 
then adjust them accordingly to ensure that each dwarf 
galaxy is detected as a single source in the segmenta¬ 
tion map, and that the full flux is being captured when 
viewing the residuals map. This segmentation map is 
also used to select the spaxels to be summed for SFR 
estimation. 

2.3. Stellar Masses 

Stellar masses are estimated according to the color- 
derived M/L ratio of West et al. (2010). The West et al. 
(2010) stellar mass estimation is a modification of the 
work done by Bell et al. (2003). The Bell et al. (2003) 
mass estimation assumes a “Diet Salpeter” IMF, how¬ 
ever the West et al. (2010) estimation modifies the Bell 
et al. (2003) estimation to a Kroupa IMF. This stellar 
mass estimation technique has been shown to be accu¬ 
rate within 20% of the actual stellar mass (Bell et al. 
2003). Photometric data for the IFU observed dwarf 
galaxies are derived from the Source Extractor measure¬ 
ments of the SDSS g, r, and i filters as detailed in Section 
2.2. With proper SDSS photometric measurements ob¬ 
tained, we utilize the West et al. (2010) stellar mass to 
light estimation: 

\og{M^/Li) = -0.222 + 0.864(^^ - r) + log(0.71) (1) 

Assuming a solar z-band absolute magnitude of 4.57 
(Sparke & Gallagher 2000), we then convert the z-band 
photometry into stellar mass estimations using the con¬ 
version outlined in West et al. (2010) along with the M/L 
ratio defined in Equation 1. 

M* = X 10-((^-4-57)/2.51) (2) 

Because these stellar mass estimations rely upon z- 
band absolute magnitudes, we require accurate distance 
estimations. For the local universe dwarf IFU observed 
galaxies, we cannot rely upon the redshift to provide ac¬ 
curate distances due to gravitational interactions within 
the local universe. Therefore, we utilize the distance es¬ 
timations from the ALFALFA catalog for all galaxies be¬ 
low a redshift of 6,000 km s~^ {z <0.02) where the errors 
on the distance estimates become comparable to those of 
Hubble flow estimations (Haynes et al. 2011). The AL¬ 
FALFA catalog distances are calculated using the local 
universe peculiar velocity model of Masters et al. (2004) 
derived from the SFI++ catalog of galaxies (Springob 
et al. 2007) and published literature values are adopted 


as appropriate. See Haynes et al. (2011) for a full ex¬ 
planation of the model used. Above a redshift of 6,000 
km s~^ we calculate the luminosity distance (D^,) from 
the SDSS redshift assuming the motions of galaxies are 
dominated by the Hubble flow. 

Using this method, stellar masses for each of the 11 
galaxies presented here are calculated to be in the range 
6.58 < log(M*) < 8.78. The full list of stellar masses 
can be found in Table 1. Uncertainties in our stellar 
masses are derived from propagating the Source Extrac¬ 
tor provided photometric uncertainties through the stan¬ 
dard Taylor series technique. 

2.4. ALFALFA/SDSS Sample Seleetion 

To compliment the IFU observed dwarf galaxy results, 
we compare this sample to galaxies selected from the 
ALFALFA HI survey. The ALFALFA team has cross- 
matched their HI detections to the SDSS database in 
order to associate the radio Hl-gas detection to optical 
photometric galaxy detections. We cross reference the 
ALFALFA a.40 database (15,855 objects) with the SDSS 
DR12 spectroscopic database and find 7,773 ALFALFA 
galaxies with SDSS spectroscopic observations. We first 
remove the 6 galaxies from the IFU observed dwarf 
galaxy sample that overlap with the ALFALFA/SDSS 
sample in order to ensure that they are not counted twice 
within our combined sample. We remove AGN contami¬ 
nation from our sample with the Kauffmann et al. (2003) 
selection criteria, leaving us with 6,151 galaxies. We re¬ 
move galaxies for which the aperture correction is ex¬ 
ceedingly large by cutting those galaxies for which the 
difference between the Petrosian r-band magnitude and 
the fiber r-band magnitude exceeds 5 mag, bringing our 
sample down to 6,137 galaxies. We also select for galax¬ 
ies with a redshift z <0.05, bringing our sample to 5,759 
galaxies. We then perform a cut to ensure that the SDSS 
spectra contain Ra emission at a S/N greater than 3 leav¬ 
ing 5,751 galaxies. After performing a final cut to exclude 
galaxies for which the Petrosian stellar mass estimations 
disagree with the MPA-JHU catalog stellar mass esti¬ 
mations by more than 0.5 dex, we have 5,436 galaxies 
in our final crossmatched ALFALFA/SDSS sample. In 
comparison, Bothwell et al. (2013) perform similar cuts, 
although they specify Ra S/N > 25 and remove galax¬ 
ies for which their two oxygen abundance estimations 
disagree, leaving them with a sample containing 4,253 
galaxies. 

2.5. Comparison Data 

To provide additional galaxies in the low stellar-mass/ 
low luminosity regime that is poorly constrained by 
the ALFALFA/SDSS sample, we also include in our 
MZR/LZR analysis data from similar surveys found in 
the literature of low stellar-mass/low luminosity galax¬ 
ies. We will utilize data from the following surveys: 

• Saintonge (2007) - Observations of 25 low Hl-mass 
ALFALFA dwarf galaxies observed using long- 
slit spectroscopy from the Double Spectrograph 
(DBSP; Oke & Gunn 1982) at Mount Palomar’s 
Hale 5m telescope. Stellar masses are derived from 
B and V band photometry using the simple stellar 
population models from Bell et al. (2003). These 
are galaxies very similar in nature to those of the 
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IFU observed dwarf galaxy sample, except that 
narrow-band Ra imaging had revealed the pres¬ 
ence of at least one bright HII region per galaxy 
making long-slit spectroscopy possible. 

• Berg et ah (2012) - A sample of low-luminosity 
galaxies selected from the Spitzer Local Volume 
Legacy survey and observed using MMT long-slit 
spectroscopy. Stellar masses are derived from B 
and K band photometry using models from Bell & 
de Jong (2001). 

• James et al. (2015) - A morphologically selected 
survey of 12 SDSS dwarf irregular galaxies, ob¬ 
served using MMT long-slit spectroscopy. Stellar 
masses are derived from Bell et al. (2003) g and r 
band photometry. 

• Survey of HI in Extremely Low-mass Dwarfs 
(SHIELD; Haurberg et al. 2015) - An ALEALEA 
derived survey of 8 galaxies observed using long- 
slit spectroscopy from the Mayall 4m telescope at 
KPNO. Stellar mass estimations come from IR (3.6 
and 4.5 jam) photometric colors using the method 
from Eskew et al. (2012). B-band magnitudes come 
from WIYN 3.5m observations. 

Due to offsets between oxygen abundance calibrations 
(Kewley & Ellison 2008), it is necessary to perform oxy¬ 
gen abundance estimations consistently in order to ob¬ 
tain meaningful comparisons between samples. There¬ 
fore, for each of the comparison samples from other stud¬ 
ies, we have recalculated oxygen abundance estimations 
(as described in Section 3.1) using the emission line fluxes 
reported within each survey. 

Similarly stellar mass estimations from all of the above 
surveys have been rescaled onto a Kroupa initial mass 
function (IME) to match the stellar mass estimations 
that we perform on both the dwarf lEU observed sample 
and the larger crossmatched ALEALEA/SDSS sample in 
Section 2.3. 

3. ANALYSIS 
3.1. Oxygen Abundance 

As a proxy for the gas-phase metallicity of a galaxy, 
we estimate the oxygen abundance using strong emis¬ 
sion line based measurements. Due to the low sur¬ 
face brightness of the lEU observed dwarf galaxy sam¬ 
ple, we are unable to detect any of the faint auroral 
lines required for direct (Tg based) oxygen abundance 
measurements. We are also limited by our wavelength 
range to using the [OHI], [Nil], and Balmer hydrogen 
emission lines to perform our oxygen abundance esti¬ 
mations. We attempted to edit the VIMOS pipeline 
source code to extend the wavelength range which would 
allow us to capture the [OH] feature required for R23 
(R23= ([OH]A3727+ [OHI]A4958,5007)/H/3) based esti¬ 
mations, however issues with interference between over¬ 
lapping pseudo-slits on the CCD proved to be a hinder- 
ance to accurate [OH] emission line measurements. 

Ideally we would use the full range of emission lines 
within our spectrum. The 03N2 based oxygen abun¬ 
dance estimation technique of Pettini & Pagel (2004) 


(hereafter PP04) is calculated using the flux of the emis¬ 
sion lines in [OHI] A5007, H/3 A4861, Ra A6563, and 
[Nil] A6583. Unfortunately, the calibration of the PP04 
03N2 based estimations becomes unreliable near or be¬ 
low 12-|-log(0/H) = 8.09 (Pettini & Pagel 2004). Eor 
completeness, we include in Appendix B a parallel anal¬ 
ysis utilizing the 03N2 oxygen abundances and find that 
our conclusions would not change. 

N2 based estimations rely only on the Ra A6563 and 
[Nil] A6583 features and are defined as 

= (3) 

N2 based oxygen abundance estimations have the ben¬ 
efit of being less sensitive to the ionization parameter 
(although not completely independent; see Morales-Luis 
et al. 2014). The PP04 N2 estimation was calibrated 
down to metallicities as low as 12-flog(0/H) = 7.48 (Pet¬ 
tini & Pagel 2004). 

Throughout this paper we will utilize the N2 oxygen 
abundance estimation as calibrated by (Denicolo et al. 
2002; hereafter D02) which was targeted towards low 
gas-phase metallicity galaxies, calibrated over a range 
7.2 < 12+log(0/H) < 9.1, and is well behaved at the 
low metallicity regime. In the D02 system the oxygen 
abundance is defined as: 

12 + log(0/H) = 9.12 + 0.73 x N2. (4) 

The blending of the emission lines due to instrumen¬ 
tal resolution in the dwarf lEU observed galaxy sample 
causes N2 based oxygen abundance estimations to be¬ 
come less accurate in low flux spaxels. However, when 
we use the integrated AoN selected spaxels for the en¬ 
tire galaxy, the integrated flux is of sufficient quality 
to reliably deblend the three Gaussian components that 
compromise [Nil] and Ra as demonstrated in Eigure 
2. Eor our dwarf galaxy lEU observed sample, we cal¬ 
culate oxygen abundance estimations in the range 7.67 
< 12+log(OH) < 8.56, indicating that our low stellar 
mass sample does indeed exhibit low gas-phase metallic- 
ity. Uncertainties in our oxygen abundance estimations 
are derived from Taylor series propagation of the un¬ 
certainties in the emission line fluxes measured in the 
lEU observed dwarf galaxy sample. The full list of oxy¬ 
gen abundances calculated for each lEU observed dwarf 
galaxy can be found in Table 1. 

3.2. Luminosity 

The B-band Luminosity-Metallicity Relation (LZR) 
has been claimed to be as strong of a gas-phase metal- 
licity indicator as the MZR when reliable distance mea¬ 
surements and appropriate photometry are used (Berg 
et al. 2012). Throughout this work, we will perform a 
parallel analysis to the EMR using luminosity in place 
of stellar mass. We calculate the B-band absolute mag¬ 
nitudes (M-b) from SDSS photometry using the relation 
from Lupton (2005). 

mB=gA 0.3130 ^{g-r)^ 0.2271. (5) 

We calculate B-band luminosities for the lEU observed 
dwarf galaxy population in the range -16.11 < M^ < 
-12.08. As detailed in Section 2.3, this requires accu¬ 
rate luminosity distances, therefore we follow the same 
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procedure and utilize ALFALFA catalog distances for 
cz <6,000 and assume a Hubble flow for cz >6,000. Un¬ 
certainties in our B-band luminosities are derived from 
propagating the Source Extractor provided photometric 
uncertainties through the standard Taylor series tech¬ 
nique. The full list of luminosities and their uncertainties 
can be found in Table 1. 

3.3. Star Formation Rates 

Star Formation Rates (SFRs) for each galaxy are de¬ 
rived from the integrated flux of Ha emission within the 
segmentation map selected spaxels of each galaxy. Spa¬ 
tial maps of the stellar continuum and Ha emission for 
each galaxy can be seen in Figure 3. The SFR for each 
galaxy is determined by first integrating the flux of every 
spaxel selected by the segmentation map selection crite¬ 
ria (the gray shaded regions in Figure 3). We use the inte¬ 
grated flux because, at the average distance of our galax¬ 
ies, each spaxel subtends 25 pc on a side and Ha based 
SFR estimations have been shown to be unreliable on 
scales less than 100-1000 pc (Kennicutt & Evans 2012), 
or when the luminosity drops below 10^^-10^^ erg s“^ 
(which corresponds to an SER of ^ 0.001-0.01 Mq yr“^). 
This would mean that several galaxies (AGC220837, 
AGC221004, AGG225852, and AGG228882) have unre¬ 
liable SER estimates on a spaxel by spaxel basis as can 
be seen in Eigure 3. 

To convert the Ha emission line measurements into 
SER estimations, we convert our Ha emission line flux 
(Ha) into luminosity (Hai,) using the luminosity dis¬ 
tance (D/,) as 

(HQ;L)obs = iTTBaDl. (6) 

where the subscript “obs” denotes the observed Ha lu¬ 
minosity. 

We need to correct for reddening to obtain the intrinsic 
Ha luminosity (HaL)int where 

= (7) 

The attenuation at a specific wavelength (Ax) is defined 

as 

Aa = k{X)E{B - V) (8) 

where E{B — V) is the broadband color excess and n{X) 

is the value of an attenuation curve at wavelength A. 
In order to obtain the reddening for each galaxy based 
on the intrinsic value of the Balmer decrement, we will 
utilize the color excess between the Ha and Hf3 emission 
lines. This can be found by substituting the wavelengths 
of Ha and H/3 into Equation 8 to find 

A{HP) - A(Ha) = n(HP)E{B -V)- n{Ha)E{B - V). 

( 9 ) 

Equation 9 could alternatively be defined as E(H/3 —Ha), 
which is analogous to E{B — V) such that 

Aim - Am = Em - m = 2.5 x log^o | | 

( 10 ) 

in which we have used the fact that the intrinsic Ha to 
H/3 ratio is 2.86 for Gase B recombination with a temper¬ 
ature T = lO^K and an electron density = 10^ cm“^ 
(Osterbrock 1989) 


By setting the right hand sides of Equations 9 & 10 
equal to each other and rearranging, we find that 


E{B - V) 


2.5 

n{Hp) - /^(Ha) 


logio 


Ram\ 
2.86 J 


( 11 ) 


Which can then be used in Equation 8 to find the color 
excess at the wavelength for Ha. We use the reddening 
curve from Galzetti et al. (2000) to obtain n{Ha) = 3.33 
and k{HP) = 4.60 for the attenuation correction: 


A(Ha) = 


2.5 X n(Ha) 

/t(H/3) - /t(HQ;) 

= 6.56 X logjo <! 


f Ha/H/3 
( 2.86 
Ram] 
2.86 J 


} 


( 12 ) 


Once we have our dust corrected Ha luminosity values, 
we apply the Hao et al. (2011) conversion 


5FR = log((HaL)int)-41.27 (13) 


We calculate star formation rates for the lEU observed 
dwarf galaxy population in the range - 2.88 < log(M 0 
yr“^)< - 1.65. The full list of star formation rates can 
be found in Table 1. Uncertainties in our star formation 
rates are propagated using standard Taylor series error 
propagation techniques from the Ha flux measurement 
uncertainties obtained via 1000 MGMG iterations. 


3.4. ALFALFA/SDSS Analysis 

In order to compare results derived for the lEU 
observed dwarf galaxy sample with the larger AL- 
EALEA/SDSS sample, we must ensure that the AL- 
EALEA/SDSS data is reduced in a manner as similar 
as possible to the lEU observed dwarf galaxy data. HI 
masses between samples are obtained from the same AL- 
EALEA database, therefore we are confident that all HI 
masses are derived similarly between the two samples. 
Eor the stellar-mass, luminosity, gas-phase metallicity, 
and SER calculations, we detail the similarities and dif¬ 
ferences between analysis of the two datasets in the fol¬ 
lowing subsections. Eigure 4 displays a comparison of 
the observables for 6 galaxies that overlap between the 
lEU observed dwarf galaxy sample and the larger AL- 
EALEA/SDSS sample. In general we find that they agree 
within 1 standard deviation, and outlying data points 
can be explained. This leads us to conclude that compar¬ 
isons between the samples are valid, however the outliers 
noted highlight the need for lEU observations and careful 
analysis in the low stellar-mass/luminosity regime. The 
sample selection criteria outlined in Section 2.4 is used to 
remove these outlier galaxies from our ALEALEA/SDSS 
sample. 


3.4.1. Stellar Mass 

Eollowing the procedure in Section 2.3, we perform the 
same calculations as in Equations 1 & 2 to obtain stellar- 
mass estimations. However, we use Petrosian magnitudes 
from the SDSS public database in place of Source Ex¬ 
tractor derived photometry. Although it would be ideal 
to perform the same Source Extractor procedure, do¬ 
ing so would be prohibitively time consuming for the 
entire ALEALEA/SDSS sample. Eor the crossmatched 
ALEALEA/SDSS sample we utilize ALEALEA catalog 
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Fig. 3. — Spatial mapping of the Ha emission line flux within a subset of the IFU fleld-of-view (Figure 1). The gray regions indicate 
the stellar continuum as selected by the segmentation maps. The blue-green pixels indicate the spaxels which pass the AoN of 5 cut. The 
intensity of the flux varies widely between galaxies, and even within galaxies the flux intensity appears to be quite patchy. Also shown for 
each galaxy is the SFR for each galaxy calculated using the flux integrated over the entire gray region showing that our galaxies cover a 
range of approximately 1.5 dex in SFR. 


TABLE 1 

Observed Properties of Dwarf Galaxies 


Galaxy 

RA Dec Distance® 

HI Mass“ 

Stellar Mass 

Luminosity 

Metallicity 

SFR 

AGG# 

hh:mm:ss.s ±hh:mm:ss Mpc (± 2.43) 

log(M0) 

log(M0) 

B Band mag 

12+log(0/H) 

log(M0 yr-l) 


AGG191702 

09:08:36.5 

+05:17:32 

8.7 

± 

2.43 

7.74 ± 0.18 

6.67 

± 

0.61 

-12.09 

± 

0.63 

7.94 

± 

0.13 

-2.76 

± 

0.34 

AGG202218 

10:28:55.8 

+09:51:47 

19.6 

± 

2.43 

7.75 ± 0.5 

8.12 

± 

0.38 

-14.95 

± 

0.28 

8.22 

± 

0.11 

-1.72 

± 

0.29 

AGG212838 

11:34:53.4 

+11:01:10 

10.3 

± 

2.43 

7.6 ± 0.19 

6.94 

± 

0.56 

-12.68 

± 

0.53 

8.05 

± 

0.14 

-2.73 

± 

0.32 

AGG220755 

12:32:47.0 

+07:47:58 

16.4 

± 

2.43 

7.18 ± 1.2 

7.76 

± 

0.43 

-13.94 

± 

0.33 

8.49 

± 

0.16 

-2.51 

± 

0.64 

AGG220837 

12:36:34.9 

+08:03:17 

16.4 

± 

2.43 

7.41 ± 0.54 

8.78 

± 

0.46 

-16.11 

± 

0.33 

8.56 

± 

0.13 

-2.2 

± 

1.28 

AGG220860 

12:38:15.5 

+06:59:40 

16.4 

± 

2.43 

7.22 ± 1.39 

7.57 

± 

0.42 

-13.92 

± 

0.33 

7.82 

± 

0.13 

-2.14 

± 

0.14 

AGG221000 

12:46:04.4 

+08:28:34 

16.5 

± 

2.43 

7.46 ± 0.83 

8.35 

± 

0.44 

-15.44 

± 

0.33 

8.35 

± 

0.05 

-1.65 

± 

0.08 

AGG221004 

12:46:15.3 

+10:12:20 

16.7 

± 

2.43 

7.66 ± 0.55 

7.98 ± 0.3 

-14.61 

± 

0.22 

8.35 

± 

0.16 

-2.23 

± 

0.22 

AGG225852 

12:59:41.9 

+10:43:40 

16.6 

± 

2.43 

7.68 ± 0.53 

7.57 

± 

0.42 

-13.89 

± 

0.33 

8.27 

± 

0.22 

-2.53 

± 

0.38 

AGG225882 

12:03:26.3 

+13:27:34 

23.6 

± 

2.43 

8.15 ± 0.3 

7.06 

± 

0.35 

-13.66 

± 

0.25 

7.95 

± 

0.18 

-2.33 

± 

0.05 

AGG227897 

12:50:04.2 

+06:50:51 

16.6 

± 

2.43 

7.44 ± 0.89 

6.58 

± 

0.45 

-12.08 

± 

0.35 

7.67 

± 

0.42 

-2.88 

± 

0.24 


Note. — “ Values obtained from the ALFALFA a.40 catalog (Haynes et al. 2011). Uncertainties in the distance are dominated by 
the local velocity dispersion measured by (Masters et al. 2005) 
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and the ALFALFA/SDSS population in order to keep 
the analysis as consistent as possible between the two 
samples. 


3.4.2. Luminosity 

Following the procedure in Section 3.2, we once again 
utilize the Petrosian magnitudes provided by the SDSS 
catalog in order to derive B-band luminosities using 
Equation 5. When we compare the results of the Pet¬ 
rosian magnitude derived luminosities to the Source Ex¬ 
tractor derived luminosities in Eigure 4, once again we 
find that AGC221004 exists considerably offset from the 
expected 1:1 correlation. As observed in Section 3.4.1, 
this is likely due to the known issues with the SDSS 
pipeline aggressively splitting up the photometry of irreg¬ 
ular HII regions. This demonstrates the need for careful 
analysis of dwarf galaxies in the low-luminosity regime. 
Errant galaxies such as AGC221004 are removed from 
the ALEALEA/SDSS sampled by the quality control cuts 
outlined in Section 2.4. 


Fig. 4. — Comparison between the measured values in galaxies 
that overlap between the IFU observed dwarf galaxy sample and 
the larger ALFALFA/SDSS sample, showing the consistency be¬ 
tween the data used to calculate observables for each of the two 
data sets. Discrepancies in the stellar-mass and luminosity values 
can be explained by the issues with SDSS Petrosian magnitude 
measurements. The poor correlation between SFR estimations can 
be explained by the patchy dwarf galaxy star formation being ob¬ 
served with a single SDSS fiber. Galaxies in blue have passed the 
selection criteria show reasonable agreement between both sam¬ 
ples, galaxies in red would be rejected from the ALFALFA/SDSS 
sample using the selection criteria outlined in Section 2.4. 

model distances for galaxies where cz < 6,000 km s“^ 
and cosmological distances otherwise. In Eigure 4 we 
find that 5 of the 6 galaxies that exist in both samples 
are within 1 standard deviation of the 1:1 correlation, 
with one significant outlier (AGG221004). 

To relieve this tension, we examine the MPA-JHU 
value added catalog (Kauffmann et al. 2003; Brinch- 
mann et al. 2004; Tremonti et al. 2004) which has stellar 
mass estimations based on spectral energy distributions 
(SED) fits to the SDSS fiber spectroscopy. These stellar- 
mass estimations utilize Bruzual & Chariot (2003) stel¬ 
lar templates and assume a Kroupa IME (Kauffmann 
et al. 2003). We find that for galaxies which overlap be¬ 
tween the lEU observed dwarf galaxy sample and the 
ALEALEA/SDSS sample, the Source Extractor based 
stellar-mass estimations agree within 7% with the MPA- 
JHU SDSS spectroscopic fiber derived total stellar mass 
estimates, suggesting that the SDSS Petrosian photom¬ 
etry is inaccurate for AGC221004. 

As outlined in West et al. (2010), roughly 25% of Pet¬ 
rosian fluxes for the SDSS population are inaccurate due 
to issues with parent-child splitting in the SDSS algo¬ 
rithms, this is consistent with what we have observed 
in Eigure 4. The MPA-JHU catalog also contains a sig¬ 
nificant number of outliers, so would not be a complete 
replacement for the SDSS petrosian magnitudes, instead 
we have implemented a quality control cut on the AL¬ 
EALEA/SDSS data in which we exclude galaxies from our 
sample for which the stellar mass estimates disagree by 
more than 0.5 dex between the Petrosian and MPA-JHU 
stellar mass estimations. With this quality control check 
in place, we continue to use the West et al. (2010) based 
mass estimations for the lEU observed dwarf galaxies 


3.4.3. Gas-phase Metallicity 

To derive the gas-phase metallicity estimations for the 
larger ALEALEA/SDSS sample, we rely upon emission 
line fluxes from SDSS single fiber spectroscopy. We apply 
the same D02 oxygen abundance estimation (Equation 
4) as was performed for the lEU observed dwarf sample. 
In Eigure 4 we find good agreement between both sam¬ 
ples. This is somewhat expected considering that the 
SDSS fibers target the brightest part of a galaxy, and 
the integrated lEU fluxes are dominated by the same 
bright areas. It is worth noting that this lends cre¬ 
dence to our emission line de-blending because the SDSS 
spectroscopy is higher spectral resolution and not signif¬ 
icantly blended. 

3.4.4. Star Formation Rate 

To determine SERs for the larger ALEALEA/SDSS 
sample, we use the Ha flux reported in the SDSS 3" 
single fiber spectroscopy database. We correct for dust 
using the same procedure as with the lEU observed dwarf 
galaxy sample. In addition, we perform an aperture cor¬ 
rection to our SER calculations. As can be seen in Eig¬ 
ure 1, the SDSS fiber often misses a significant amount 
of flux in the emission regions of nearby galaxies. There¬ 
fore, to more accurately compare our lEU observations 
which cover the entire star forming region to the SDSS 
observations, we apply the aperture correction technique 
outlined in Hopkins et al. (2003) 

—0-4(rpetro—'^fiber) (U) 

where Upetro denotes the SDSS r-band Petrosian flux, and 
Ffiber denotes the SDSS r-band flux within the fiber. 

As observed in Eigure 4, the calculated star formation 
rates may vary by as much as 0.5 dex in SER between the 
two samples, emphasizing the need for lEU observations 
to collect all of the Ha flux. This is to be expected 
considering the patchy star forming nature of the lEU 
observed dwarf galaxy sample. The patchy nature of 
star formations also means that this aperture correction 
is imperfect, as seen is Eigure 4. These facts indicate 
that ALEALEA/SDSS SERs are likely to be unreliable 
for galaxies requiring a large aperture correction. This 
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Fig. 5.— The mass-metallicity relation (top) and the luminosity- 
metallicity relation (bottom) for the IFU observed dwarf galaxy 
sample, along with several low stellar mass galaxy observations and 
the ALFALFA/SDSS sample. The scatter in both the MZR (LZR) 
fit around the bin means is 0.05 dex (0.03 dex). For reference we 
plot a sun symbol at solar metallicity, showing that all of the IFU 
observed dwarf galaxy oxygen abundances are below solar. The 
yellow regions indicate the calibration limits of the D02 N2 oxygen 
abundance estimation. We have used the published emission line 
fluxes from other authors to estimate the oxygen abundance in the 
same D02 system. The MZR appears to hold down to the lowest 
stellar masses observed (log(M>,c) = 7.25). 

necessitates the sample selection cut where we exclude 
galaxies for which the difference between the Petrosian 
r-band magnitude and the fiber r-band magnitude differ 
by more than 5 mag as described in Section 2.4. 


4. RESULTS 

4.1. Mass-Metallicity Relation and 
Luminosity-Metal licity Relation 

The mass-metallicity relation (MZR) can be seen in 
Figure 5 along with homogenized observations from ear¬ 
lier dwarf/low luminosity galaxy surveys. Similar to the 
process outlined in Mannucci et al. (2010), we bin the 
ALFALFA/SDSS crossmatched sample by stellar mass, 


into bins 0.25 dex wide. We then fit a 4 degree polyno¬ 
mial to the mean values of the bins and find the following 
relationship: 

12 + logiOjH) = 8.756 + 0.219 x m - 0.119 x 
—0.052 X — 0.003 x 

(15) 

where m = log(M^)-10. The MZR is fit to data from the 
literature collected samples, the ALFALFA/SDSS cross- 
matched sample, and the IFU observed dwarf galaxy 
sample. In order to ensure robust results, we require 
that there are at least 21 galaxies per stellar mass bin. 
We use the limit of 21 galaxies per bin in order to be 
consistent with the analysis performed in Bothwell et al. 
(2013). In total there are 15 stellar-mass bins covering 
the stellar-mass range 7.25 < log(M*) < 11.0. The resid¬ 
ual Icr scatter of every galaxy to the MZR fit in Equation 
15 is 0.11 dex, similar to the 0.1 dex scatter reported for 
the MZR of Tremonti et al. (2004). The Icr scatter of the 
mean values in each stellar-mass bin is 0.05 dex. The Icr 
scatter of the means is the metric that we will use to 
compare scatter between various metallicity relations. 

The LZR can also be seen in Figure 5 along with ho¬ 
mogenized observations from similar dwarf/low luminos¬ 
ity galaxy surveys. We bin the combined literature, AL¬ 
FALFA/SDSS, and IFU observed sample into bins 0.52 
dex wide. We then fit a 4 degree polynomial to these 
mean values and find the following relationship: 

12 + log{0/H) = -18.6 - 7.23 x Mb - 0.724 x Mb^ 

-0.032 X Mb^ - 0.0005 x Mb'‘. 

(16) 

The LZR is fit to 15 luminosity bins over a range —13.7 < 
M^) < — 21.0 in order to ensure that there are at least 
21 galaxies per luminosity bin. The residual Icr scatter 
of every galaxy to the LZR fit in Equation 16 is 0.14 
dex. The residual Icr scatter of the mean values to the 
LZR in each luminosity bin is 0.03 dex. Again, this is 
the metric that we will use to compare scatter between 
various metallicity relations. 

4.2. MZR/LZR with SFR and HI Mass Dependance 

We test whether or not a fundamental metallicity rela¬ 
tion explains the enhanced scatter of the MZR and LZR 
by following a similar procedure as Mannucci et al. (2010) 
and Bothwell et al. (2013). We bin the crossmatched AL- 
EALEA/SDSS data into SER bins of 0.85 log(MQ yr“^) 
and HI mass bins of 0.85 log(M 0 ). We use larger HI 
mass and SER bins than Bothwell et al. (2013) to en¬ 
sure that we have at least 21 objects per bin in our lower 
stellar mass/luminosity regimes. Eigures 6 & 7 show the 
result of binning the MZR and the LZR in terms of both 
SER and HI mass. A clear separation exists between 
the HI mass bins shown in the top of Eigures 6 & 7, as 
well as the SER bins in the bottom of Eigures 6 & 7. 
The numerical values of the means used to produce the 
solid lines in Eigures 6 & 7 can be found in Appendix C. 
Using these larger HI mass and SER bins, we are able 
to extend the ALEALEA/SDSS means relation down to 
stellar mass bins as low as log(M*) = 7.75. 

Also shown in Eigures 6 & 7 are the linear least squares 
fits to the lEU observed dwarf galaxies. The dwarf galaxy 
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Fig. 6.— The mass-metallicity relation as seen in Figure 5 color- 
coded by HI mass (top) and SFR (bottom). The stars indicate 
individual observations within our IFU observed sample of dwarf 
galaxies. The dashed lines indicate linear least squares fits to these 
points, and the shaded regions indicate the Icr standard deviations 
to these fits. The solid curves indicate the mean values of the AL¬ 
FALFA/SDSS sample, separated into HI mass bins (top) and SFR 
bins (bottom). Shown in the background are the ALFALFA/SDSS 
points which are binned to produce the color-coded means. We 
find little overlap between the HI mass bins (top) and SFR bins 
(bottom). 

sample of this study is not large enough to produce sim¬ 
ilarly binned mean relationships for dwarf IFU observed 
data alone, therefore we instead perform least squares 
optimized linear fits to the dwarf IFU data to demon¬ 
strate that it is consistent with the FMR. The numerical 
values for the linear fits and their la standard deviations 
can be seen in Table 2. 

The IFU observed dwarf galaxy sample overlaps with 
the larger ALFALFA/SDSS sample by ^ 1 dex, provid¬ 
ing a valuable consistency check. The added benefit of 
the IFU observed dwarf galaxy sample is that it is able 
to test the fundamental metallicity relation down to stel¬ 
lar masses 2 orders of magnitude lower than observed in 
Mannucci et ah (2010) and Bothwell et ah (2013). 

4.3. Fundamental Metallicity Relations 


Fig. 7. — The luminosity-metallicity relation as seen in Figure 
5 color-coded by HI mass (top) and SFR (bottom). The stars 
indicate individual observations within our IFU observed sample 
of dwarf galaxies. The solid curves indicate the mean values of 
the ALFALFA/SDSS sample, separated into HI mass bins (top) 
and SFR bins (bottom). Shown in the background are the AL¬ 
FALFA/SDSS points which are binned to produce the color-coded 
means. We find little overlap between the HI mass bins (top) and 
SFR bins (bottom). 

In order to test whether the IFU observed dwarf galaxy 
sample is constant with the offsets in SFR and HI Mass 
observed for the larger ALFALFA/SDSS crossmatched 
sample, we examine the fundamental metallicity relation¬ 
ships. We also introduce nomenclature for the luminosity 
equivalent of the FMR as the fundamental metallicity lu¬ 
minosity (FML) relation, which can be dependent upon 
SFR (FMLsfr) or HI-ps mass (FMLri). The FMR 
and FML are a projection of three-dimensional param¬ 
eter space into two dimensions to minimize the scatter 
observed. 

4.3.1. MZR SFR Dependence (FMRsfr) 

Mannucci et al. (2010) define a variable which com¬ 
bines SFR with stellar mass such that 

/i, = log(M*)-Mog(SFR). (17) 

They find that a = 0.32 minimizes the scatter in the 
FMRsfr relation, whereas Bothwell et al. (2013) find 


log(Hl Mass) [log(M©)] SFR (M© yr 
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Fig. 8.— Left: la scatter in the means of the FMRsfr as a function of a. The optimal scatter found by Mannucci et al. (2010) is 
indicated by the dashed green vertical line, and the optimal alpha value found by Bothwell et al. (2013) is indicated by the dashed red 
vertical line. We find the lowest scatter in the FMRsfr when a = 0.28 Right: Icr scatter in the means of the FMRhi as a function of /3. 
The dashed red vertical line indicates the optimal /3 value from Bothwell et al. (2013) (No ^ minimum exists from Mannucci et al. (2010) 
as they did not study the FMRhi )• We find the lowest scatter in the FMRhi when ^ = 0.33. 


TABLE 2 

Linear Fits to Dwarf IFU SFR & HI Mass Bins 


SFR Bin 

Stellar Mass 
Slope 

Y-Intercept 

Icr Scatter 

-2.4 to -1.55 

0.44 ± 0.01 

4.7 ± 0.7 

0.15 

-3.25 to -2.4 

0.5 ± 0.0 

4.59 ± 0.22 

0.13 

HI Mass Bin 

Slope 

Y-Intercept 

Icr Scatter 

6.75 to 7.6 

0.3 ± 0.01 

5.85 ± 0.85 

0.22 

7.6 to 8.45 

0.24 ± 0.0 

6.34 ± 0.26 

0.1 


Luminosity 



SFR Bin 

Slope 

Y-Intercept 

Icr Scatter 

-2.4 to -1.55 

-0.27 ± 0.0 

4.15 ± 0.64 

0.13 

-3.25 to -2.4 

-0.28 ± 0.0 

4.54 ± 0.5 

0.16 

HI Mass Bin 

Slope 

Y-Intercept 

Icr Scatter 

6.75 to 7.6 

-0.15 ± 0.0 

6.07 ± 0.85 

0.24 

7.6 to 8.45 

-0.13 ± 0.0 

6.27 ± 0.32 

0.11 


Note. — Linear fit results for the dotted lines as 
seen in Figures 6 & 7. The Icr scatter is the width of the 
shaded regions in Figures 6 & 7. 


a = 0.28 is the optimal value. 

We perform the same tests on our larger crossmatched 
ALFALFA/SDSS sample in order to find the value of 
a minimizes the FMRsfr scatter. To find the optimal 
value of d, we project the ALFALFA/SDSS data onto 
the /i-metallicity plane and bin our data using the same 
SFR bins used for Figure 6. Each SFR binned dataset 
is binned a 2nd time in /i space similar to the stellar- 
mass binning performed in Section 4.2 and we once again 
perform the same cut of 21 galaxies per /r-SFR bin. We 
take the mean metallicity of each /r-SFR bin and fit a 4th 
degree polynomial to the binned mean metallicity values. 
This process is repeated for values of a ranging from -1 
to 2. 

The dataset used for this procedure includes only 
the larger ALFALFA/SDSS dataset. We tested using 


the combined ALFALFA/SDSS and IFU observed dwarf 
data and found the same scatter per alpha value tested. 
This is due to the fact that there are only 11 data points 
in the IFU observed data set being added to the consid¬ 
erably larger ALFALFA/SDSS sample. Also producing 
the FMRsfr using only the larger ALFALFA/SDSS data 
means that the test of whether or not the IFU observed 
dwarf galaxy data fits the same FMRsfr is completely 
independent. 

If Q/ = 0 provides the optimal value, that would indi¬ 
cate that the MZR has the lowest scatter obtainable and 
that taking into account the SFR is unnecessary. Also 
of note, if a = 1 provides the optimal value, that would 
indicate that the specific star formation rate (M*/SFR) 
would be the best projection to use. The scatter obtained 
for each value of a can be seen in the left hand side of 
Figure 8. We find that a value of a = 0.28 provides 
the lowest scatter for the FMRsfr (Table 3). Using this 
a value we then reproduce the FMRsfr plot using our 
optimal projections for Figure 9. 

For the FMRsfr relation seen in Figure 9, we find that 
although the 2nd lowest (light blue) SFR bin overlaps 
between the IFU observed dwarf galaxy sample and the 
larger crossmatched ALFALFA/SDSS sample, the low¬ 
est (dark blue) SFR bin is not consistent between the 
IFU observed dwarf galaxy sample and the crossmatched 
ALFALFA/SDSS sample. In fact the lowest SFR bin is 
offset Icr above the FMRsfr- This is suggestive of a 
breakdown of the FMRsfr for the lowest SFR galaxies. 
It is possible that our SFRs are slightly inaccurate in 
this range considering they are near the limit discussed 
in Kennicutt & Evans (2012). As star formation rates 
approach this region, they should be treated as upper 
limits of the true SER. 

4.3.2. MZR HI-Mass Dependence (FMRhi) 

Bothwell et al. (2013) defines an equation similar to 
Equation 17 for the EMRhi relation: 

rip = log(M*) - l3{\og{MHi) - 9.80). (18) 
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Fig. 9.— FMRhi (top) and FMRsfr (bottom) as calculated 
using the a (top) and ^ (bottom) values found to minimize scatter 
(Figure 8). We find that the FMRhi relation is consistent down 
to the lowest Hl-mass bin, where as the lowest SFR bin in the 
FMRsfh is offset Icr above the FMRsfr- 


Wherein they find that a value of (3 = 0.35 minimizes 
the scatter in the gas-phase metallicity-?] plane. 

To find the optimal value of /3, we project the AL¬ 
FALFA/SDSS data onto the 77 -metalhcity plane and bin 
our data using the same HI mass bins used for Figure 
6 . Each HI mass binned dataset is binned a 2nd time 
in 77 space similar to the stellar-mass binning performed 
in Section 4.2 and we once again perform the same cut 
of 21 galaxies per 77 -HI mass bin. We take the mean 
metallicity of each 77 -HI mass bin and fit a 4th degree 
polynomial to the binned mean metallicity values. This 
process is repeated for values of p ranging from -1 to 2 . 

If /3 = 0 provides the optimal value, that would indi¬ 
cate that taking into account the HI mass is unnecessary. 
Also of note, if /3 = 1 provides the optimal value, then 
the ratio of Hl-gas mass to stellar mass would be the best 
projection to use. The scatter obtained for each value of 
P can be seen in the right hand side of Figure 8 . We 
find that a value of p = 0.32 provides the lowest scatter 
for the FMRhi (Table 3). Using this P value we then 
produce the FMRhi plot using our optimal projections 
for Figure 9. 


For the FMRhi relation, we find that our linear fits to 
the IFU observed dwarf galaxy sample are, within the 
Icr uncertainties, consistent with the binned means of 
the larger ALFALFA/SDSS crossmatched survey. This 
suggests that the FMRhi relation is likely to continue 
down to stellar masses as low as 10 ^’^ M^. 

4.3.3. LZR SFR Dependence (FMLsfr) 

We define a fundamental luminosity-metallicity rela¬ 
tion (FMLsfr) using the star formation rate as: 

ec = (Mb) - Clog(SFR) (19) 

To find the optimal value of we project the AL¬ 
FALFA/SDSS data onto the e-metallicity plane and bin 
our data using the same SFR bins used for Figure 7. 
Each SER binned dataset is binned a 2nd time in e space 
similar to the luminosity binning performed in Section 
4.2 and we once again perform the same cut of 21 galax¬ 
ies per e-SER bin. We take the mean metallicity of each 
e-SFR bin and fit a 4th degree polynomial to the binned 
mean metallicity values. This process is repeated for val¬ 
ues of C ranging from -2 to 1 . 

The scatter obtained for each value of can be seen in 
the right hand side of Figure 10. We find that a value of C 
= -0.35 provides the lowest scatter for the FMLsfr (Ta¬ 
ble 3). Using this value we then produce the FMLsfr 
plot using our optimal projections for Figure 11. 

For the FMLsfr relation seen in Figure 11, we find 
that although the 2nd lowest (light blue) SFR bin over¬ 
laps between the IFU observed dwarf galaxy sample and 
the larger crossmatched ALFALFA/SDSS sample, the 
lowest (dark blue) SFR bin is not consistent between 
the IFU observed dwarf galaxy sample and the cross- 
matched ALFALFA/SDSS sample. In fact the lowest 
SFR bin is offset Icr above the FMLsfr derived from the 
ALFALFA/SDSS sample. This is suggestive of a break¬ 
down of the FMLsfr relation for the lowest SFR galax¬ 
ies. However, as discussed in Section 4.3.1, it is possible 
that our SFRs are inaccurate in this range considering 
they are near the limit discussed in Kennicutt & Evans 
( 2012 ). 

4.3.4. LZR HFMass Dependence (FMLhi) 

We define a fundamental luminosity-metallicity rela¬ 
tion (FMLhi) using the star formation rate as: : 

Ct = (Mh) — 7(log(MH-7) — 9.80). (20) 

To find the optimal value of 7 , we project the AL¬ 
FALFA/SDSS data onto the ^-metallicity plane and bin 
our data using the same HI mass bins used for Figure 
7. Each HI mass binned dataset is binned a 2nd time 
in ^ space similar to the luminosity binning performed in 
Section 4.2 and we once again perform the same cut of 21 
galaxies per ^-HI mass bin. We take the mean metallic- 
ity of each ^-HI mass bin and fit a 4th degree polynomial 
to the binned mean metallicity values. This process is 
repeated for values of 7 ranging from -2 to 1 . 

The scatter obtained for each value of 7 can be seen in 
the right hand side of Figure 8 . We find that a value of 7 
= -0.50 provides the lowest scatter for the FMLhi (Table 
3). Using this 7 value we then produce the FMLhi plot 
using our optimal projections for Figure 11. 

For the FMLhi relation, we find that our linear fits to 
the IFU observed dwarf galaxy sample is, within the Icr 
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Fig. 10.— Left: la scatter in the FMLsfr relation as a function of (. We find the lowest scatter in the FMLsfr relation when = —0.35 
Right: Icr scatter in the means of the FMLhi as a function of 7 . We find the lowest scatter in the FMLhi relation when 7 = -0.50. 
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Fig. 11 .— FMLhi (top) and FMLsfr (bottom) as calculated 
using the ^ and 7 values found to minimize scatter (Figure 10). 
We find that the FMLhi relation is consistent down to the lowest 
Hl-mass bin, where as the lowest SFR bin in the FMLsfr is offset 
la above the FMLsfr- 


TABLE 3 

Best Fitting Fundamental Relation 
Parameters 


Ma 

a 

log(M,) - alog(SFR) 

0.28 

np 

/3 

log(M.)-/3(log(Mjf/)-9.80) 

0.32 


c 

(Mb) - Clog(SFR) 

-0.35 


7 

{Mb) — 7 (log(MH/) — 9.80) 

-0.50 


Note. — The values that produce the 
lowest scatter in Figures 8 & 10. 


uncertainties, consistent with the binned means of the 
larger ALFALFA/SDSS crossmatched survey, suggesting 
that the FMLhi relation is likely to continue down to 
luminosities as low as = —12. 

5. DISCUSSION 

As can be seen in Figure 5, the MZR holds across 5 
orders of magnitude in stellar mass, with a la scatter of 
0.05 dex. There is an upturn in the MZR around M* 
= 10^ Mq which can be explained in the context of the 
FMR. Galaxies with stellar masses below 10^ Mq are 
likely to have low star formation rates and low Hl-gas 
masses which correlates with higher than average metal- 
licity for a specific stellar mass. For a galaxy with 10^ 
Mq, for it to exist in the highest HI mass bin, it would 
have to have an Hl-gas mass 2 orders of magnitude larger 
than the stellar-mass. 

The LZR also holds over 7 magnitudes, with a la scat¬ 
ter of 0.03 dex. A similar upturn can be seen in the LZR 
around M^ = -15. Again this can be explained by the 
expected low Hl-mass and/or low SFR of low luminosity 
galaxies. The lower scatter in the LZR supports the hy¬ 
pothesis of Berg et al. (2012) that accurately measured 
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luminosities perform equally as well as stellar-masses. 

To provide insight into the most physically motivated 
components of a fundamental metallicity relation, we ex¬ 
amine the Icr scatter of the means in the larger AL- 
FALFA/SDSS data set for each permutation of mass, 
luminosity, SFR, and HI mass used to form a possible 
fundamental relation. Using the scatter of the means 
around each quartic fit as the evaluation criteria, we find 
that the lowest scatter of the means is in the combination 
of stellar mass and HI mass binning (FMRhi) for which 
we obtain a mean scatter of 0.01 dex, the other three 
relationships (FMRsfr, FMLhi, FMLsfr) have mean 
scatters of 0.02 dex. This suggests that the FMRhi may 
be the most physically significant fundamental relation. 

5.1. FMRhi 

In examining the FMRhi relation, seen in Figure 9 
we find that the dwarf IFU sample is consistent within 
1 standard deviation with the FMRhi relation observed 
within the larger ALFALFA/SDSS sample. This suggests 
that the FMRhi continues to hold down to stellar masses 
as low as M* = 10^-^ Mq. Our results suggest that the 
physical mechanism responsible for the FMRhi must be 
active across the entire stellar mass range from 10^’^ to 
10^0 M0 

Bothwell et al. (2013) reach the conclusion that the 
FMRhi is more physically motivated due to the reduced 
scatter and physical motivation of the FMRhi- They 
also noticed a lack of saturation in the high mass end of 
the FMRhi- In out analysis, which uses slightly different 
stellar-mass and Hl-mass bins, we do observe saturation 
in stellar mass bins above log(M*) = 10.0 in Figure 6. We 
find the FMRhi relation to be more physically motivated 
due to its behavior on the low mass end (log(M*) < 9.0) 
in which the IFU observed dwarf galaxy data is consistent 
with the larger ALFALFA/SDSS dataset. The la scatter 
of the FMRhi binned mean metallicities being the lowest 
of the 4 permutations supports this hypothesis. 

Inflows of pristine gas diluting the metal content would 
explain why the FMR has an HI mass dependance. Kop- 
pen & Edmunds (1999) showed that the gas-phase metal- 
licity of a galaxy can indeed be reduced if the gas infall 
rate is larger than the rate at which gas is converted 
into stars. Our results suggest that inflows of pristine 
gas continue to drive the FMRhi down to stellar masses 
as low as 10^-^ Mq. This agrees with recent models by 
Grpnnow et al. (2015) which show that pristine gas rich 
mergers are partially responsible for the scatter in the 
FMR. 


5.2. FMR SFR 

Unlike the FMRhi, the FMRsfr (Figure 9) does not 
appear to hold across our entire sample. The lowest SFR 
bin in the IFU observed dwarf galaxy sample is offset la 
higher than the FMRsfr as determined using the larger 
ALFALFA/SDSS sample. If our results are accurate near 
the calibration limits of the SFR and oxygen abundance, 
they suggest a breakdown in the physical mechanism re¬ 
sponsible for the FMRsfr ia very low SFR populations. 

Possible explanations for the FMRsfr have focused 
on outflows of metal-enriched gas caused by the intense 
winds of young star-forming regions within a galaxy. 
Simulations by Mac Low & Ferrara (1999) have shown 
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Fig. 12.— FMRsfr with the two lowest SFR bins combined into 
one in order to compensate for the possible underestimation of SFR 
in the lowest SFR bin of Figure 9. We see that in this scenario, 
the FMRsfr is consistent across the stellar mass range observed. 

that supernova winds are capable of expelling metals ef¬ 
ficiently enough to contribute to the MZR. We find that 
this is likely to be true down to an SFR of -2.4 log(M^ 
yr“^), however, below that it is possible that stellar 
winds are not sufficiently strong enough to remove metal- 
enriched gas, and galaxies have a higher gas-phase metal- 
licity than would be expected by the FMRsfr- Although 
limitations in our measurements could also explain this 
apparent breakdown. 


5.3. FMLhi 

We find that the FMRhi reiation, seen in Figure 11, 
is consistent within la between the IFU observed dwarf 
gaiaxy popuiation and the iarger ALFALFA/SDSS pop- 
uiation. These resuits are similar to those observed 
for FMRhi however we find that the scatter of the 
binned ^- 0.5 mean metaiiicities around the quartic fits is 
iarger than the anaiogous fits for FMRhi suggesting that 
steiiar-mass is the more physicaiiy motivated parameter. 

5.4. FMLsfr 

The FMLsfr reiation is inconsistent in both the high 
and iow regions, as can be seen in Figure 9. Sim- 
iiar to the FMRsfr the iowest SFR bin is > Icr offset 
from the FMLsfr relation as determined by the larger 
ALFALFA/SDSS sample. The highest SFR bin also devi¬ 
ates by as much as Icr, suggesting that the high mass SFR 
bin is also poorly fit by the best fitting 77 in the FMLsfr- 
The inconsistencies observed within the FMLsfr reiation 
provide further evidence that Hl-gas mass may be the 
more physicaiiy motivated component in a fundamental 
metallicity relation. 

5.5. Ha Line Flux Limitations 

In order to explain the apparent breakdown in the 
FMRsfr and FMLsfr relations for the lowest SFR bins, 
we consider the possibility that our Her line flux limit may 
be causing inaccurate SFR estimations. In comparing 
Her flux measurements between the AoN selected binned 
spectra and the segmentation map selected binned spec¬ 
tra, we find the segmentation map selection captures ap- 
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proximately 33% more Ha flux indicative of low SFR re¬ 
gions being missed. It is also possible that metallicities 
in the fainter, low SFR, galaxies are over estimated due 
to limitations in measuring low Nil A6583 fluxes, placing 
them preferentially above the FMRsfr (See Appendix 
A). 

To test the hypothesis that our lowest SFR bins are 
underestimated, and our lowest oxygen abundances are 
overestimated, we re-examine the FMRsfr considering 
the two lowest SFR bins as one bin. Figure 12 demon¬ 
strates that the combined bin is consistent with the 
FMRsfr- Therefore, it is possible that the FMRsfr 
is indeed valid over the same stellar mass range as the 
FMRhi however, care must be taken when obtaining IFU 
observations to ensure that the target depth is sufficient 
to accurately estimate the SFR and oxygen abundance. 
To determine whether line flux measurement limitations 
or physical limitations are the cause of the la offset of 
the lowest SFR bin, we would require further observa¬ 
tions with longer integration times. 


6 . CONCLUSION 

In order to extend the fundamental relations down to 
stellar masses as low as 10^-^ and luminosities as low 
as -12 mag, we utilize a sample of IFU observed dwarf 
galaxies in combination with a larger ALFALFA/SDSS 
crossmatched population of galaxies. Whenever possible, 
we perform the same analysis on both populations in 
order to ensure that results are consistent between the 
two data sets. 

When comparing galaxies that overlap between the 
SDSS and IFU observed samples, we find that IFU ob¬ 
servations are necessary to accurately measure the star 
formation rates of the dwarf irregular galaxies presented 
here due to the patchy nature of their star forming re¬ 
gions. However, gas-phase metallicites appear to be ad¬ 
equately measured using a single SDSS fiber (Figure 4). 

Using a sample of galaxies selected from the AL¬ 
FALFA blind HI survey, in combination with galaxies 
observed with the VIMOS IFU spectrograph and a se¬ 
lection of long-slit observed galaxies, we investigate the 
mass-metallicity relation (MZR) with a particular focus 
on the low stellar mass/low luminosity regime. We find 
that the MZR continues down to stellar masses as low as 
107.25 (Figure 5). We find that the MZR has a la 
scatter in the means around a quartic fit of 0.05 dex. 

We find using the ALFALFA/SDSS crossmatched pop¬ 
ulation that the fundamental metallicity relation as a 
function of SFR (FMRsfr) exhibits its lowest scatter in 
the means (la = 0.02 dex) when a = 0.28. We also find 
using the IFU observed dwarf galaxy population that 
the FMRsfr relations appears to break down in the low¬ 
est SFR bins (Figure 9). However limitations of strong- 
line metallicity estimations, and Ha based SFRs in low 
stellar-mass dwarf galaxies may provide an alternative 
explanation for this apparent breakdown (Section 5.5). 

We find using the ALFALFA/SDSS crossmatched sam¬ 
ple that the fundamental metallicity relation (FMRhi) 
exhibits the lowest scatter in the means (la = 0.01 dex) 
when (3 = 0.32. We also find that the FMRhi relations 
are both consistent between the IFU dwarf population 


and the ALFALFA/SDSS population across the entire 
Hl-gas mass range. 

We find that the luminosity-metallicity relation (LZR) 
continues down to B-band luminosities as low as -14 mag. 
We find that the LZR has a la scatter in the means 
around the quartic fit of 0.03. The lowest scatter in the 
FMLhi it is 0.02 dex (for C = -0.35). Similarly we find 
that the lowest scatter in the FMLsfr is 0.02 dex for 
7=-0.50. 

When comparing our sample of dwarf galaxies to the 
Mannucci et al. (2010) and Bothwell et al. (2013) analy¬ 
sis, we find that our results are consistent with theirs in 
that our derived a and jS values agree ±0.04. In agree¬ 
ment with the conclusions of Bothwell et al. (2013), we 
find that the FMRhi is more physically motivated than 
FMRsfr based on the most significant reduction of scat¬ 
ter. 

In summary, the FMRhi appears to be the most phys¬ 
ically significant driver of the fundamental metallicity 
relation, suggesting that inflows of pristine gas would be 
a possible explanation for the fundamental relations ob¬ 
served. 
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APPENDIX 

APPENDIX A - ABILITY TO RECOVER METALLICITY 

Due to the spectral resolution and instrumental dispersion of the VIMOS LR Blue grism, some emission lines are 
severely blended together, as can be seen if Figure 2. Our science results depend on the [Nil] and Ha emission lines 
which are separated by ^20 A. To test our ability to recover the emission line fluxes of the blended [Nil] and Ha 
emission lines within our spectra, we create a simulated galaxy with known emission line ratios and test our ability to 
recover them using our Python-based Gaussian fitting routines. 

To create the simulated galaxy, we first take a pipeline reduced data cube from one of the galaxies in our sample 
that was too faint to be detected by the IFU spectrograph (AGC220261) and find a region of the field of view which 
contains only sky with little instrumental artifacts to interfere. We also check to ensure that none of the galaxy that 
was originally intended to be observed is located near this region. We then generate a simulated galaxy and place it 
within this region. This simulated galaxy is then run through our full data reduction pipeline to measure the recovered 
emission line ratios. This process is repeated several times to test for metallicity dependance of our ability to recover 
emission line ratios. 

To ensure that we generate physically realistic simulated galaxies, we draw our input emission line ratios from the 
Berg et al. (2012) survey of low-luminosity galaxies. For every galaxy in the Berg et al. (2012) sample, we use the 
long slit measured emission line ratios to create a simulated galaxy with the same emission line ratios. We scale the 

brightest central spaxel in the simulated galaxy to have an H/3 emission line flux of 3.0 x 10“^^ erg cm“^ s“^ A , 
which corresponds to the average H(3 flux in the brightest spaxel of each galaxy in the dwarf IFU sample. The flux of 
the simulated galaxy decreases exponentially with radius. Each emission line is given a dispersion equal to the average 
instrumental dispersion of 7.5 A. 
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Fig. 13.— Example of the Ha emission observed in a simulated galaxy as recovered by our Python-based Gaussian fitting routines. Also 
shown is the recovered Ha flux after integrating the spectra of each spaxel shown. 


An example of the Ha emission map for one of our simulated galaxies can be seen in Figure 13. Emission lines for a 
sample fiber of our simulated galaxy can be seen in Eigure 14 which can be compared to Eigure 2 visually to confirm 
that the simulated galaxies being created are physically realistic. 
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Fig. 14. — Left: Input simulated spaxel for our simulated galaxy. Blue curves indicate the simulated spectral data, red lines indicates 
the linear continuum fit done in the first stage of our Python-based Gaussian fitting routines, and green points indicate the fitting window 
that was used to determine the linear fit. The fitting window excludes residuals from skylines and the large internal reflection feature at 
~ 6300A. Right: Three Gaussian fits of the Ha and [Nil] emission lines as found by our routines. The routines are able to successfully 
resolve the 3 blended emission lines. 























Recovered Metallicity Recovered Metallicity 


19 


APPENDIX B - COMPARISON BETWEEN N2 AND 03N2 BASED CALIBRATIONS 
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Eig. 15.— Showing the ability of our Python-based Gaussian fitting routines to recover metallicities in various systems after creating a 
simulated galaxy. In the top row, the x-axis shows the input gas-phase metallicity based on the input line flux ratios, and the y-axis shows 
the result measured using our routines. In the bottom row, the x-axis shows the direct metallicity that was used to create the simulated 
galaxy, and the y-axis shows the gas-phase metallicity measured using our routines. We find that the D02 oxygen abundance calibration 
is more accurate than the PP04 based systems, however PP04 03N2 is more precise than the D02 based systems. 


We simulate every galaxy in the Berg et al. (2012) sample using the method outlined in Appendix A to test both 
our ability to recover input strong-line emission ratios, and to compare our recovered metallicity estimations to direct 
metallicity measurements. The top row of Figure 15 shows the comparison between input and recovered metallicities 
in each of the three calibrations: PP04 03N2, PP04 N2, and D02 N2. We find that our Python-based Gaussian 
fitting routines are able to successfully recover line fluxes within 1 standard deviation down to the lowest metallicity 
simulated. The larger error bars in N2 line ratios of the low-metallicity simulated galaxies measured are due to the 
amplitude of the [Nil] lines being comparable to the noise of the spectra. 

The bottom row of Figure 15 shows the gas-phase metallicity measured directly along with the recovered gas-phase 
metallicity. Although the 03N2 based estimation has smaller error bars due to the increased information from the 
OIII] and H/3 emission line information, it is less accurate than the N2 based estimations due to the saturation of the 
OIIIJ/H/3 ratio known to occur in 03N2 based oxygen abundance estimations at low gas-phase metallicity (Pettini & 
Page! 2004). 

We find based on the bottom row of Figure 15 that D02 N2 oxygen abundance estimations are better correlated with 
direct metallicity measurements than PP04 03N2 estimations for low-metallicity galaxies such as those presented in 
this study. We find that PP04 03N2 provides more precise, but less accurate observations, whereas D02 N2 provides 
more accurate, but less precise observations for the lowest metallicity galaxies. 

For the sake of completeness, we have reproduced all of the figures from the main text which contain oxygen 
abundance estimations using the 03N2 based oxygen abundance estimation from Pettini & Pagel (2004): 


(Bl) 


12 + log(0/H) = 8.73 - 0.32 x 03N2 
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where 


03N2 = log 


[OIII]X5007/RP\ 
[NII]X65S3/Ra J 


(B2) 


We find that in general, the trends observed in the paper and the conclusions of the paper continue to hold, especially 
on the high-metallicity end, however the oxygen abundance floor of the 03N2 based estimations cause the conclusions 
to be less solid than as observed in the N2 based results of the main text. 
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Fig. 16. — Same as Figure 5 however using the PP04 03N2 oxygen abundance calibration in place of the D02 N2 calibration. Left: 
The mass-metallicity relation for the IFU observed dwarf galaxy sample, along with several low stellar mass galaxy observations and the 
ALFALFA/SDSS sample. Right: The luminosity-metallicity relation for the IFU observed dwarf galaxy sample, along with several low 
stellar mass galaxy observations and the ALFALFA/SDSS sample. The scatter in both the MZR and the LZR fit around the bin means is 
0.05 dex. For reference we plot a sun symbol at solar metallicity, showing that all of the IFU observed dwarf galaxy oxygen abundances are 
below solar. The yellow regions indicate the calibration limits of the D02 N2 oxygen abundance estimation. We have used the published 
emission line fluxes from other authors to estimate the oxygen abundance in the same PP04 03N2 system. The mass-metallicity relation 
appears to flatten out in the low stellar-mass/low luminosity regime. 
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Fig. 17.— Same as Figure 6 however using the PP04 03N2 oxygen abundance calibration in place of the D02 N2 calibration. The 
mass-metallicity relation as seen in Figure 16 color-coded by HI mass (left) and SFR (right). The stars indicate individual observations 
within our IFU observed sample of dwarf galaxies. The dashed lines indicate linear least squares fits to these points, and the shaded regions 
indicate the Icr standard deviations to these fits. The solid curves indicate the mean values of the ALFALFA/SDSS sample, separated into 
HI mass bins. Shown in the background are the ALFALFA/SDSS points which are binned to produce the color-coded means. We And little 
overlap between the SFR bins (right). 
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Fig. 18.— Same as Figure 7 however using the PP04 03N2 oxygen abundance calibration in place of the D02 N2 calibration. The 
luminosity-metallicity relation as seen in Figure 16 color-coded by HI mass (left) and SFR (right). The stars indicate individual observations 
within our IFU observed sample of dwarf galaxies. The solid curves indicate the mean values of the ALFALFA/SDSS sample, separated 
into HI mass bins. Shown in the background are the ALFALFA/SDSS points which are binned to produce the color-coded means. We find 
little overlap between the SFR bins (right). 
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Fig. 19.— Same as Figure 9 however using the PP04 03N2 oxygen abundance calibration in place of the D02 N2 calibration. FMRhi 
(left) and FMRsfr (right) as calculated using the a and /3 values found to minimize scatter (Figure 8). We find that the FMRhi relation 
is consistent down to the lowest HI-mass bin, where as the lowest SFR bin in the FMRsfr is offset Icr above the ALFALFA/SDSS derived 
FMRsfr- 
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Fig. 20.— Same as Figure 11 however using the PP04 03N2 oxygen abundance calibration in place of the D02 N2 calibration. FMLhi 
(left) and FMLsfr (right) as calculated using the C and 7 values found to minimize scatter (Figure 10). We find that the FMLhi relation 
is consistent down to the lowest Hl-mass bin, where as the lowest SFR bin in the FMLsfr is offset Icr above the FMLsfr- 
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APPENDIX C - SFR AND HI MASS MEAN BINS 

To correspond to the linear fits, we provide a set of tables showing the mean metallicity values in each one of the 
SFR and HI mass binned data for both the MZR and LZR. These values could be used to reproduce the solid curves 
of the SFR or HI-mass bins in Figures 6 & 7. 
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TABLE 4 

ALFALFA/SDSS SFR Binned Mean Metallicity Values 


Stellar Mass 

12+log(0/H) 

12+log(0/H) 

12-blog(0/H) 

12-blog(0/H) 

12-blog(0/H) 

Bin Range 

Unbinned SFR Range 

SFR -2.6 to -1.7 

SFR-1.7 to -0.8 

SFR -0.8 to 0.1 

SFR 0.1 to 1.0 

7.75 to 8.0 

8.22 

8.22 




8.0 to 8.25 

8.23 

8.25 

8.17 

- 

- 

8.25 to 8.5 

8.29 

8.36 

8.26 

- 

- 

8.5 to 8.75 

8.35 

8.41 

8.33 

- 

- 

8.75 to 9.0 

8.42 

8.52 

8.43 

8.28 

- 

9.0 to 9.25 

8.5 

- 

8.53 

8.42 

- 

9.25 to 9.5 

8.58 

- 

8.62 

8.55 

- 

9.5 to 9.75 

8.67 

- 

8.70 

8.67 

8.58 

9.75 to 10.0 

8.74 

- 

8.79 

8.74 

8.70 

10.0 to 10.25 

8.78 

- 

- 

8.79 

8.77 

10.25 to 10.5 

8.81 

- 

- 

8.83 

8.80 

10.5 to 10.75 

8.83 

- 

- 

8.84 

8.82 

10.75 to 11.0 

8.82 

- 

- 

- 

8.82 

Stellar Mass 

12+log(0/H) 

12-blog(0/H) 

12-blog(0/H) 

12-blog(0/H) 

12-blog(0/H) 

Bin Range 

Unbinned HI Mass Range 

HI Mass 7.6 to 8.45 

HI Mass 8.45 to 9.3 

HI Mass 9.3 to 10.15 

HI Mass 10.15 to 11.0 

7.75 to 8.0 

8.22 

8.30 




8.0 to 8.25 

8.23 

8.22 

8.25 

- 

- 

8.25 to 8.5 

8.29 

8.37 

8.27 

- 

- 

8.5 to 8.75 

8.35 

- 

8.35 

8.35 

- 

8.75 to 9.0 

8.42 

- 

8.44 

8.40 

- 

9.0 to 9.25 

8.5 

- 

8.53 

8.48 

- 

9.25 to 9.5 

8.58 

- 

8.63 

8.57 

- 

9.5 to 9.75 

8.67 

- 

8.72 

8.67 

- 

9.75 to 10.0 

8.74 

- 

8.78 

8.74 

8.67 

10.0 to 10.25 

8.78 

- 

8.78 

8.78 

8.75 

10.25 to 10.5 

8.81 

- 

- 

8.81 

8.80 

10.5 to 10.75 

8.83 

- 

- 

8.83 

8.82 

10.75 to 11.0 

8.82 

- 

- 

8.82 

8.82 

Luminosity 

12+log(0/H) 

12-blog(0/H) 

12-tlog(0/H) 

12-blog(0/H) 

12-blog(0/H) 

Bin Range 

Unbinned SFR Range 

SFR -2.6 to -1.7 

SFR-1.7 to -0.8 

SFR -0.8 to 0.1 

SFR 0.1 to 1.0 

-15.26 to -14.74 

8.29 

8.27 




-15.78 to -15.26 

8.33 

8.34 

8.28 

- 

- 

-16.3 to -15.78 

8.37 

8.40 

8.33 

- 

- 

-16.81 to -16.3 

8.38 

8.43 

8.36 

- 

- 

-17.33 to -16.81 

8.47 

8.56 

8.46 

8.43 

- 

-17.85 to -17.33 

8.52 

8.52 

8.53 

8.50 

- 

-18.37 to -17.85 

8.59 

- 

8.61 

8.57 

- 

-18.89 to -18.37 

8.66 

- 

8.69 

8.65 

8.71 

-19.41 to -18.89 

8.73 

- 

8.81 

8.73 

8.72 

-19.93 to -19.41 

8.77 

- 

- 

8.78 

8.77 

-20.44 to -19.93 

8.81 

- 

- 

8.82 

8.80 

-20.96 to -20.44 

8.81 

- 

- 

- 

8.80 

Luminosity 

12+log(0/H) 

12-blog(0/H) 

12-tlog(0/H) 

12-blog(0/H) 

12-blog(0/H) 

Bin Range 

Unbinned HI Mass Range 

HI Mass 7.6 to 8.45 

HI Mass 8.45 to 9.3 

HI Mass 9.3 to 10.15 

HI Mass 10.15 to 11.0 

-14.74 to -14.22 

8.24 

8.24 




-15.26 to -14.74 

8.29 

8.30 

8.20 

- 

- 

-15.78 to -15.26 

8.33 

8.32 

8.31 

- 

- 

-16.3 to -15.78 

8.37 

- 

8.35 

- 

- 

-16.81 to -16.3 

8.38 

- 

8.37 

8.37 

- 

-17.33 to -16.81 

8.47 

- 

8.48 

8.46 

- 

-17.85 to -17.33 

8.52 

- 

8.56 

8.51 

- 

-18.37 to -17.85 

8.59 

- 

8.61 

8.59 

- 

-18.89 to -18.37 

8.66 

- 

8.70 

8.66 

8.64 

-19.41 to -18.89 

8.73 

- 

8.79 

8.73 

8.72 

-19.93 to -19.41 

8.77 

- 

- 

8.77 

8.77 

-20.44 to -19.93 

8.81 

- 

- 

8.80 

8.82 

-20.96 to -20.44 

8.81 

“ 

“ 

8.81 

8.80 


Note. — Binned means used to produce the solid lines observed in Figures 6 & 7. 
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APPENDIX D - TABLE OE ELUX VALUES 

We provide our emission line measurements for the integrated spectra from each galaxy’s IFU observations so that 
others may use them for their own analysis. 


TABLE 5 

Observed Emission Line Eluxes of Dwarf Galaxies 


Galaxy AGC# 4861 A [OIII] 5007 A Ha 6563 A [Nil] 6583 A 


AGG191702 

AGG202218 

AGG212838 

AGG220755 

AGG220837 

AGG220860 

AGG221000 

AGG221004 

AGG225852 

AGG225882 

AGG227897 


82.1 ± 11.3 

56.3 ± 8.3 
47.0 ± 6.3 

10.7 ± 3.7 

12.3 ± 3.3 

60.7 ± 5.1 

125.1 ± 7.2 

38.1 ± 13.5 

21.8 ± 9.4 

52.4 ± 7.5 
14.3 ± 0.6 


247.4 ± 25.5 

223.4 ± 28.0 
174.3 ± 21.3 

23.4 ± 7.7 

35.1 ± 9.8 

321.2 ± 23.3 

230.6 ± 12.5 

69.9 ± 23.4 
48.8 ± 19.1 

111.2 ± 12.8 

74.2 ± 1.3 


249.6 ± 31.8 

223.7 ± 31.1 

151.1 ± 20.5 
43.0 ± 13.9 
65.1 ± 15.9 

212.2 ± 14.9 
492.5 ± 29.3 

126.2 ± 43.0 

66.5 ± 28.3 

142.0 ± 21.4 

41.7 ± 1.0 


6.0 


2.3 

13.1 


4.0 

5.2 


2.2 

5.9 


2.3 

11.2 


3.9 

3.6 


1.5 

43.5 


6.3 

11.0 


4.1 

4.6 


2.6 

3.6 


1.9 

0.4 


0.6 


Note. — All line fluxes measurements are in units of 
10 “^® erg cm“^ s“^ A . These are the line fluxes used to calculate 
oxygen abundance estimations throughout this work. 
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